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In supergravity models with low supersymmetry breaking scale the gravitinos can be superlight, 
with mass in the 10 -6 eV to few keV range. In such a case, gravitino emission provides a new 
cooling mechanism for protoneutron stars and therefore can provide constraints on the mass of 
a superlight gravitino. This happens because the coupling to matter of superlight gravitinos is 
dominated by its goldstino component, whose coupling to matter is inversely proportional to the 
scale of supersymmetry breaking and inceases as the gravitino mass decreases. Present observations 
therefore provide lower limits on the gravitino mass. Using the recently revised goldstino couplings, 
■ we find that the two dominant processes in supernova cooling are e + e~ — > GG and 7 + e~ — ■» e~GG. 

They lead to a lower limit on the supersymmetry breaking scale As from 160 to 500 GeV for core 
temperatures 30 to 60 MeV and electron chemical potentials 200 to 300 MeV. The corresponding 
lower limits on the gravitino mass are .6 — 6 x 10 -6 eV. 
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I. INTRODUCTION 

In supergravity models with a low supersymmetry breaking scale, the gravitino (G) mass is in the superlight range 
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of 10 cV to few keV's since it is given by the formula mg ~ 2^3M — (where As is the scale of supersymmetry 
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breaking and Mpi is the Planck mass that characterises the gravitational interactions.) Any information on the 
gravitino mass therefore translates into knowledge of one of the most fundamental parameters of particle physics, the 
. scale at which supersymmetry breaks. 

There is an advantages in studying the gravitino's properties when its mass is in the superlight range because 
it can be easily emitted in astrophysical processes such as supernova cooling ffl, neutron star cooling etc. This 
adds new cooling mechanisms to the already known ones for supernovae, i.e. the usual neutrino emission, which 
the observation of neutrinos from SN1987A seems to have confirmed. Any additional process can therefore afford a 
maximum luminosity of roughly 10 52 ergs/sec. This will then lead to constraints on the parameters that describe 
the coupling of the gravitinos to matter in the supernovae. Furthermore since for light gravitinos one has G M ~ 

i\J\^rd fl Xi the superlight gravitino coupling is dominated by the coupling of the goldstino to matter, which is 

inversely proportional to the supersymmetry breaking scale-squared F = A|. For values of As in the 100 GeV to TeV 
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range, the goldstino coupling strengths to matter are of the same order as the ordinary weak interactions. The gravitino 
emission in astrophysical settings such as the supernovae can therefore be competitive with the neutrino emission 
process. Observed supernova neutrino luminosity by the 1MB and Kamiokande groups [|| and its understanding in 
terms of the standard model of the supernova Q therefore allows us to set lower limits on As and hence on m^. 

The supernova and other astrophysical constraints on the gravitino mass were first studied in two recent papers 
. The first paper considered the class of models where the superlight gravitino is the only superlight particle in 
the model with its superpartners having masses in the GeV range whereas the last two papers considered the smaller 
subclass of models |(| where the gravitino is also accompanied by superlight scalar and pseudoscalar particles. In this 
paper we will focus on the first class of models. 

In Ref. (I), gravitino couplings to matter suggested in Ref. Q] were used. These couplings have recently been 
criticized in two papers fjjjjioll and a new set of matter couplings have been proposed for the class of supersymmetry 
models where the scalar and pseudoscalar partners of the gravitino are heavier (e.g. in the multi GeV range) than 
the gravitino. Since the temperature dependence of gravitino emission rates are very different if one uses the new 
set of Feynman rules for gravitino matter couplings, it is necessary to revisit these bounds again. It is the goal of 
this paper to use the revised Feynman rules for matter gravitino coupling to calculate gravitino emission rates in 
supernovae, obtain the bounds on the supersymmetry breaking scale, and from them the lower limit on the gravitino 
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mass[] Our result is that the basic processes that dominate the energy loss via gravitino emission are e + e — > GG 
(called annihilation process below) and 7 + e~ — > e~GG, (called Compton process below) whereas the process found 
to dominate in Ref. M i.e. 77 — > GG is found to make a negligible contribution. The physically interesting lower 

limit on the to^ remains ~ 3 x 1CP 6 x ( lor^GeV ) ^ e ^ (where M is a model dependent parameter, which can lie 
anywhere from 50 to 250 GeV). This bound is qualitatively somewhat better than the one found in Ref. §. 

This paper is organized as follows: in section II, we note the various possible processes that can contribute to energy 
loss from supernovae via G emission and, using simple dimensional analysis, obtain the temperature dependence of 
the emissivity for the different processes and give qualitative arguments to isolate the processes that dominate the 
emissivity; we calculate the emissivity (Q) for the Compton process in sec. Ill and for the annihilation processes in 
sec. IV and derive a lower bound on Ag and vtiq from supernova; in section V we discuss the energy loss from neutron 
stars. 



II. DOMINANT MECHANISMS FOR ENERGY LOSS VIA GRAVITINO EMISSION 

Since gravitinos are superpartners of gravitons, which have universal coupling to matter, we expect them to couple 
to all matter in pairs. Specifically in the case when the gravitino is superlight, the dominant couplings arise from the 
coupling of its longitudinal mode jll). Goldstinos, like Goldstone bosons, must be derivatively coupled. Combining 
this property with the constraints of supersymmetry, Luty and Ponton Jlpf have derived the form of the effective 
coupling between matter and goldstinos. The couplings of interest to us are given by |l0| 

^IXX = a 4 3^X1 X-f^f (1) 

^ = ^d^(D^) (2) 

C = Z ^d»xi u *AF flvA (3) 
A s 

where we have denoted the goldstino by \ an d the generic gaugino by A,4 for the gauge field A. In the 7XX coupling 
the mass parameter M is highly model dependent and sensitive to the details of the supersymmetry breaking sector 
as well as the nature unification grou P0. In Ref. p| a value of 43 GeV is suggested assuming the typical slepton mass 
to be 200 GeV and the squark mass to be 400 GeV. But it is perfectly possible even in the context of gauge mediated 
models to have very different parameters e.g. slepton masses of 200 GeV and squark masses to be a TeV jl^]. We 
have therefore kept this as a free parameter and express all our results in terms of this number. 

We are now ready to discuss the detailed processes that can contribute to the energy loss in in a supernova. There 
are the following classes of processes: (i) ones that involve both nonrelativistic particles (e.g. neutrons and protons) 
and relativistic particles such as e^, 7 and \\ (ii) that involve only relativistic particles; and finally, (iii) the ones that 
involve the decay of the plasmon i.e. 7* — > GG. 

A typical process of type (i) is the analog of the modified URCA process for neutrino emission with neutrinos 
replaced by the gravitinos : NN — > NN\X- Even though the gravitino emission graph involves the exchange of 
a virtual photon, the derivative couplings present in the gravitino vertex imply that the photon-gravitino-gravitino 
coupling in momentum space goes like ki ■ = (fc 2 — 2m 2 = ,)/2 where k = k\ + &2 and ki are the gravitino momenta. 

For superlight gravitinos, the effective virtual photon induced e + e~ — > GG amplitude is like a pointlike four Fermi 
coupling with Gf/V% replaced by e 2 M 2 /8As- One can therefore use this feature to calculate the gravitino pair 



x As this paper was being prepared for publication, a paper by Grifols et al [g] appeared which used the revised Feynman 
rules to obtain a lower limit on the gravitino mass. While our final results are very similar, we find two processes dominating 
gravitino emission process whereas Ref. || considers only one of them. 

2 For instance in grand unified theories with simple groups, at the GUT scale M 2 = 0. Therefore at low energies, this will 
give additional suppression factors 
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emissivity from the known the URCA processes O]. One can conclude that the emissivity Q, defined as the amount 
of energy emitted from the supernova per unit volume per unit time (E 5 in natural units) is given by 

Qnn * (1647r 3 /4725)a 2 a 2(M/A2)VT 8 (4) 

Let us now consider several examples of the class (ii) processes. They are (a) e + e~ — * GG; (b) 7 + e~ — > e~ + G + G; 
(c) 7 + 7 — > GG. Let us first focus on the process (c) which is the weakest of the three. Since this process 
involves exchange of a photino, its coupling strength is given by ~ e 2 /(A|Af\ ). Therefore from simple dimensional 
arguments we conclude that the emissivity for this process will go like ~ (a 2 /AgM| )T 15 . Let us now apply similar 
dimensional arguments to the process (b). In this case the Feynman diagrams that contribute are given in Fig. 1. 
Note that all exchange particles are also light. Therefore the only dimensional parameters are provided by the photon- 
gravitino coupling strength M /A|. Using dimensional analysis, one gets the following temperature dependence for Q: 

~ a 3 (jx^j T 9 . Here, we have ignored the effects of Pauli exclusion principle in this estimate; this will be addressed 

below. But assuming that this effect does not lead to a drastic change in temperature dependence, we see that it 
clearly dominates over the process (c) since supernova core temperatures are of order 30 to 60 MeV. 

Process (a) is more subtle since the electrons are degenerate in the supernova core. As a result of this the positron 
density is suppressed compared to that of the electron by by (T//i) 3 e~ M / T , where fi is the chemical potential for the 
electrons. Therefore even though the naive dimensional arguments imply that for this process, the emissivity goes 

like a 2 (^j^ T 9 and is therefore more dominant than the process (b), we will see that in actual practice they are 
comparable. 

The plasmon decay has been considered in || and is found to be negligible. We therefore do not discuss it here. 



III. CALCULATION OF THE EMISSIVITY VIA THE COMPTON PROCESS 



The generic formula for the emissivity Q for a process with initial particles and momenta denoted by i a and pi a 
and final state particles and momenta denoted by f a and pf a is: 

q= / n ^ Sr2 3 ff (2 " )4 ' 3 "'" 3 " f (5) 

5\Y, Pla - Ep /a )II ia /(p iB )n /a (l - f{pu))\Mi t f?(E n + E f2 ) 

In the equation above, Mif denotes the matrix element for the process responsible for gravitino emission; Efij 2 are 
the energies of the gravitinos; f(p) are the thermal distributions for the various particles involved in the process. For 
example the thermal distribution of the electrons and positrons are given by: 

/ eT = l/(e^ + 1) (6) 

where /Lt e is the chemical potential for the electrons. Similarly for the photons we have fj(p) = l/(e9 — 1). For the 
gravitinos, we choose the /'s to be zero since they do not get a chance to thermalize in the supernova. 

Let us apply this formula to the Compton process which is given by the two Feynman diagrams in Fig.l. Let us 
denote the initial and final momenta as follows: 

7 (<z) + e(p 1 )^e(p 2 ) + G(fc 1 ) + G(fc 2 ) (7) 

In evaluating |M| 2 for this process, we note that the 7GG vertex is common to both the direct and the cross 
diagrams. Therefore, we evaluate that absolute square first as follows: 

ii , -/^^'(*-h-*>»lW < 8 » 



where 



y a aa= e ^^(Wv(k 2 ) (9) 
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Using its gauge invariance, one can write this as 

\V\ 2 = A{k a k - g^k 2 ) (10) 
where A is easily evaluated to be 4p 6 4 ^{ . The matrix element Mj / is given by 

Mif = e V%S(p 2 )[7 Q ( 7 • {q+pt) - to)- V + 7^(7 • (pa - <z) - m)~ V]«(Pi) (H) 
One can then combine the above equations to write the final formula for emissivity as follows: 

Qcom P = (||) 4 / A ^f 3g ^/(^)/-r(g)(l - /(P2))* a [E?*i/«<] (12) 

where at = (2q ■ pi) 2 , a.2 = (2q ■ P2) 2 and 03 = — 2g • piq ■ P2 and 

Xi = 32(2to 4 + m 2 (2q ■ pt - q -p 2 -p% ■ P2) + q -p%q ■ P2) (13) 
X 2 = 32(2to 4 + m 2 (q ■ p 1 - 2q ■ p 2 - P\ ■ p 2 ) + q • Pig ■ P2) 
X 3 = 16m 2 (q ■ p 2 - q • Pi) + 32pi ■ p 2 (q ■ Pi - q ■ P2 - Pi ■ P2 + 2m 2 ) 

This integral was evaluated by Monte carlo method to obtain the emissivity as a function of the supernova core 
temperature and the parameters M and Ag. Multiplying by the volume of the supernova (with radius 10 kilometers), 
we obtained the luminosity which was then set less than 10 52 ergs/sec. The bound on the parameter Ag for M = 100 
GeV are given table I. 

We also wish to point out that we have checked plasma screening effects on our result by redoing the calculation 
with the propagator k 2 replaced by k 2 + fc 2 ; (with k p i as given in Braaten and Segel [Q in the relativistic limit). We 
find that these effects are well below 5% level as might be expected from the ki ■ ^2 factor in Eq. (9). 



IV. CONTRIBUTION OF ELECTRON POSITRON ANNIHILATION TO EMISSIVITY 

Using similar methods as in section III, we have calculated the emissivity in this case and find 

Qann = 8a 2 (M / K 2 fT A e~ ^b{p / T) / lb** (14) 

where b(y) = (5/6)e y y~~ 5 (F^ + F^F^) where F^(y) — J dxx" 1 ^ 1 / (1 + e x±y ). Our expression agrees with that 
calculated in the Ref. ||. We have adopted their notation in Eq. (14). We have evaluated the integrals in the 
above expression numerically and have obtained lower bounds on the As using the requirement that the luminosity 
in gravitinos have an upper bound of 10 52 ergs/sec. The results for this case are comparable to the limits derived 
from the Compton cooling case and both cases are collected in Table I. We have varied the core temperature from 
30 to 60 MeV and used two typical values for the electron chemical potential of 200 and 300 MeV. Again, since the 
dependence of the luminosity on the supersymmetry breaking scale goes like the eighth power of temperature, for 
each case the bigger of the two numbers is the actual bound. We see that the best bound obtains for the case when 
the core temperature is assumed to be 60 MeV as expected and the chemical potential is 200 MeV and is 

A s > 500 GeV (15) 

For the case where T c — 50 MeV and p = 300 MeV, the bound is As > 300 GeV for both the Compton and the 
annihilation cases. Combining these two together we get the bound to be As > 2 1 / 8 x 300 = 325 GeV for this choice 
of supernova parameters. Our bound is actually slightly weaker than the bound given in || using only the e + e~ 
annihilation mechanism. Note that we have a higher value for M. 

Finally, as was first pointed out in pj, the above considerations assume that the supernova is transparent to the 
gravitinos once they are emitted. To check this, we have to make sure that the mean free path of the grvitino, once 
emitted is longer than the radius of the supernova (i.e. 10 km). In Ref. §, the mean free path was calculated 
using photon-gravitino scattering. Now of course due to the revised Feynman rules, the whole picture is different. 
The main contribution to opacity comes from gravitino scattering off protons and electrons. The mean free path for 
A s = 300 GeV due to proton scattering has been calculated in || and we do not repeat their calculation. They find 
that for As < 220 GeV, the gravitinos get trapped. Using the techniques employed in [fL5[ , they have shown that the 
emissivity remains too large until As < 70 GeV, so that the excluded range of As is between 70 GeV and 300- 500 
GeV depending on the choice of core temperature and the electron chemical potential. This can be translated to a 
forbidden range for the gravitino mass of 10~ 7 eV < < .6 — 6 x 10~ 6 eV. This bound is very similar to the bound 
obtained in 
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V. NEUTRON STAR COOLING CONSTRAINTS 



Neutron star cooling is an attractive area in which to test non-standard particle physics models, in principle, because 
the standard model of neutron cooling, based on the "modified URCA process," gives insufficient cooling to match 
current observations. For a recent review and refernces to the literature sec |l6[ . Briefly, for neutron stars of ages 
between about 10 2 and 10 5 years, the standard model (slow cooling) is one of neutrino emission from an isothermal 
superfluid core with an energy gap in the 100 keV range. However data from two pulsars of ages around 10 4 years 
show temperatures a factor of 5 or so below the standard model prediction. Time is measured by P/(2^£) where 
P is the pulsar period while temperature is measured from the (X-ray) black body spectrum (keeping in mind that 
the surface temperature is believed to be roughly the square root of the interior temperature (in eV)). For earlier 
times only a few upper bounds are available so that there is no good evidence on the temperature dependence of the 
mechanism that is providing the extra (fast) cooling. However, considerable extra cooling is required since the URCA 
cooling rate is proportional to the eighth power of the temperature. 

Various mechanisms to provide the extra cooling have been proposed. These include a smaller superfluid energy 
gap, meson condensates, and nucleon dissociation into quarks under the high central pressure. Were gravitino emission 
the source of the extra cooling, that fact would provide a measurment of gravitino coupling rather than a bound. 
Thus we turn to the question of whether that possibility can be entertained. 

As noted earlier, the effective coupling of the gravitino pair to electrons and quarks is via the photon exchange - 
however for supcrlight gravitinos, it reduces effectively to a four Fermi interaction like the weak Fermi coupling except 
that the coupling only occurs with electrically charged particles. One channel for neutron star cooling that involves 
gravitino emission is therefore via a neutral current like coupling with the coupling constant given by q a e 2 M 2 / (4A|). 
Thus the relevant process is the analog of the neutral current URCA process n + p — >n+p + GG. The temperature 
dependence of this process will be T 8 as in the URCA process. So it can compete with the neutrino cooling of neutron 
stars in the appropriate range for the c oup ling parameters. Assuming that neutrino luminosity describes the cooling 
of neutron stars long after their birth |lq| , we assume that luminosity via gravitinos is of same order of magnitude. 
(We choose this arbitrarily in view of the scant data on neutron stars at the moment.) 

As is well known, in the case of the neutron stars the neutral current driven URCA process is about a factor of 
30 lower |l3|,|l^] than the corresponding charged current process, of which a factor of four comes from the coupling 
and the rest from phase space. Since our couplings are given, we demand that gravitino cooling be less than the full 
neutrino URCA process luminosity. This leads us to the bound that 

4naM 2 „ . r 

~^Af- ~ GfI ^ (16) 

This implies that As > 200 GeV, which is comparable to the bounds obtained from the supernova. Considerably 
larger coupling, contrary to the supernova bounds derived above, would be required for gravitino emission to provide 
an explanation for the fast neutron star cooling observed. 



VI. CONCLUSION 



In conclusion, we have revisited the issue of supernova and neutron star constraints on the supersymmetry breaking 
scale and the ensuing constraint on the mass of the superlight gravitino using the revised Feynman rules for the 
goldstino coupling to matter in a large class of supersymmetric models with a low SUSY breaking scale. We find that 
the lower limits on the supersymmetry breaking scale As are in the range of 200 to 500 GeV and the resulting lower 
limit on the gravitino mass is .6 — 6 x 10~ 6 eV. One must however remember that while the form of the gravitino 
coupling to matter is universal, there is a model dependent parameter M that characterises the strength of the 
coupling. In any case, these bounds are comparable to the collider bounds obtained by various authors flS[ ]. As has 
been noted earlier [|J , the corresponding constraints on the gravitino mass are much more severe in models where the 
gravitino is also accompanied by superlight scalar/or pseudoscalar particles. 
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Table I 



T c 




Ag nm from Qann 


A rnm from Q Comp 


30 


200 


165 


175 


30 


300 


135 


160 


40 


200 


250 


230 


40 


300 


225 


240 


50 


200 


335 


300 


50 


300 


305 


300 


60 


200 


420 


500 


60 


300 


400 


385 



Table caption: The lower bounds on the supersymmetry breaking scale As in GeV derived from the two dominant 
processes in the supernova for various values of core temperature T c in MeV and electron chemical potential \i e in 
MeV. We have chosen M = 100 GeV. 

Figure caption Feynman diagram for Compton cooling where the lower solid line denotes the electron, the wavy 
lines denote the photon and upper external solid lines denote the gravitino pair. 
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